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ABSTRACT 

Using complete samples of steep-spectrum quasars, we present evidence 
for a correlation between radio and optical luminosity which is not caused 
by selection effects, nor caused by an orientation dependence (such as 
relativistic beaming), nor a byproduct of cosmic evolution. We argue that 
this rules out models of jet formation in which there are no parameters in 
common with the production of the optical continuum. This is arguably 
the most direct evidence to date for a close link between accretion onto a 
black hole and the fuelling of relativistic jets. The correlation also provides 
a natural explanation for the presence of aligned optical/radio structures 
in only the most radio luminous high-redshift galaxies. 



1 INTRODUCTION 

It is now established that quasars can be divided 
into two physically distinct classes: radio-loud quasars 
(RLQs) and radio-quiet quasars (RQQs) (Peacock, 
Miller & Longair 1986; Kellermann et al. 1989; Miller, 
Peacock & Mead 1991; Miller, Rawlings & Saunders 
1993; Wilson & Colbert 1995). Relativistic jets are 
a universal feature of RLQs (e.g. Bridle et al. 1994) 
and are also probably associated with at least some 
RQQs (e.g. Miller et al. 1993); the difference between 
the classes is thus not whether the quasars can form 
relativistic jets, but related instead to the fraction of 
the total power output channelled along them in a 
bulk kinetic form (see also Rawlings 1994). 

It is now also established that at least some ra- 
dio galaxies harbour obscured quasar nuclei (e.g. An- 
tonucci 1993; Antonucci, Hurt & Kinney 1994; Dey & 
Spinrad 1996, Ogle et al. 1997). However, the popular 
notion that the probability of obscuration is a strong 
function of the angle between the jet axis and the 
line-of-sight (e.g. Barthel 1989; Antonucci 1993) due 
to the ostensible presence of a dusty molecular torus, 
is still debated, and may apply only to a restricted 
range of radio luminosity and/or redshift (see e.g. 



Lawrence 1991; Jackson & Rawlings 1997 and refs. 
therein). Accepting this orientation-based unification 
scheme, however, it is possible to combine RLQs and 
at least some radio galaxies into a single 'radio-loud' 
category, and estimate the power in the photoionising 
ultraviolet continuum Q p hot (which is hidden in the 
case of radio galaxies) from the luminosity of the nar- 
row emission lines. By doing this in the 3CRR sample 
(Laing, Riley & Longair 1983) Rawlings & Saunders 
(1991) inferred that 'radio-loud' objects have bulk 
powers Qbuik in their jets of the same order as Qphot, 
whereas in the case of RQQs (and their possible ob- 
scured counterparts; see Lonsdale, Smith & Lonsdale 
1989) the fraction of the power output channeled into 
jets is 10 3 times lower (Miller et al. 1993; Rawl- 
ings 1994). The strongly differing fractions seem to 
be a fundamental difference between radio-loud and 
radio-quiet active galactic nuclei. 

Despite this difference there is some evidence 
that within each class of quasar there is some intrinsic 
connection between Qbuik and Qphot- If one accepts 
that radio luminosity (L ra d) and [O ill] line luminos- 
ity (Lro nil) are crude indicators of these variables 
then the separate, but roughly parallel, loci of the 
radio-loud and radio-quiet objects in the L ra d versus 
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L[o in] plane (e.g. Rawlings 1994) hints at a correla- 
tion between Qbuik and Qphot whether the jet is in 
a radio- loud quasar (with Qbuik ~ Qphot) or a radio- 
quiet quasar (with Qbuik "C Qphot)- This gives rise to 
the hope that these relations might constrain models 
for the fuelling of relativistic jets in all types of active 
galaxies (e.g. Rawlings & Saunders 1991; Falcke et al. 
1995). 

However, this is far from being a universally ac- 
cepted view. In both radio-loud and radio-quiet ob- 
jects, the correlation between Qbuik and the inferred 
Qphot can be influenced, or even caused by selection 
effects in samples used to date, as can the apparently 
similar magnitudes of Qbuik and Q p hot in radio-loud 
objects. We will discuss this in more detail in section 
3.4. 

There are further objections, on physical rather 
than statistical grounds. For example, Dunlop & Pea- 
cock (1993) have suggested that I/ ra d-i[o in] corre- 
lations in radio-loud objects are more plausibly ex- 
plained by models in which both luminosities are en- 
hanced through interactions between radio jets (and 
lobes) and a dense environment (see also Baum, Zirbel 
& O'Dea 1995 and refs. therein). There are certainly 
some cases of both radio-loud objects (e.g. Lacy & 
Rawlings 1994) and radio-quiet objects (e.g. Axon et 
al. 1989) where at least some of L[o in] may be at- 
tributable to power supplied by the jet. A key ques- 
tion, therefore, is whether one can find a direct link 
between Qbuik and Qphot rather than inferring the 
latter indirectly from L[o mj. 

The optical continuua and nuclear emission lines 
of RLQs are the direct probes of Qphot required. Dun- 
lop & Peacock (1993) cite the lack of a correlation be- 
tween the optical continuum and radio luminosities of 
RLQs as an argument in favour of their environmental 
interpretation of the I/ ra d-i[o in] correlation. In fact 
the existing evidence on the radio-optical correlation 
for RLQs is mixed (Peacock, Miller & Longair 1986; 
Browne & Murphy 1987; Neff, Hutchings & Gower 
1989; Miller, Peacock & Mead 1990; Miller et al. 1993 
- we will review these contradictory claims in section 
3.4) and, in practice, confused by a number of issues. 

Firstly, Peacock, Miller & Longair (1986) showed 
that one must not mix RLQs and RQQs since do- 
ing so can produce an apparently universal (but spu- 
rious) radio-optical correlation. The authors argued 
that RLQs have a minimum optical luminosity to ex- 
plain why 3C quasars are among the brightest op- 
tically, and why the RLQs in the optically-selected 
Bright Quasar Survey (Schmidt & Green 1983) are 
among the most luminous radio sources. This alone 
would yield a spurious correlation when comparing 
RLQs with RQQs over a wide optical range. (In fact 
our new data rules out this scenario, as we discuss in 
section 3.4.) 

Secondly, the strong orientation dependence of 
the optical and radio continuua in flat-spectrum core- 
dominated RLQs (e.g. Jackson et al. 1989; Baker 



1996), can lead to a radio-optical correlation. This 
follows because radio and optical luminosities can be 
simultaneously enhanced in these objects by relativis- 
tic beaming of synchrotron radiation. 

Thirdly, although there are hints of cou- 
pled optical-X-ray-radio processes in ostensibly un- 
beamed RLQs (Browne & Wright 1985; Browne & 
Murphy 1987), the samples used to date either lack 
the spectroscopic redshifts which would distinguish 
radio luminosity dependence from evolution, or have 
serious and unquantifiable selection biases. 

In this paper we present the results of a new 
study of the radio-optical correlation for quasars in 
which we have attempted to implement lessons learnt 
from the studies referenced above. We confined our 
attention to one of the distinct classes of quasars, the 
RLQs, since we believe their radio luminosities are a 
more straightforward indicator of Qbuik (which, with 
suitable radio data, it is possible to estimate in in- 
dividual quasars; Rawlings & Saunders 1991; Rawl- 
ings 1993). A recent study of the radio-optical cor- 
relation in RQQs can be found in Lonsdale et al. 
1995. To reduce the effects of relativistic beaming 
and optical continuum anisotropy we have also cho- 
sen to concentrate on steep-spectrum RLQs (here- 
after SSQs), to the exclusion of flat-spectrum, core- 
dominated quasars. Section 2 outlines the selection 
and study of our new complete sample of SSQs. To 
distinguish between L ra d-dependence and evolution 
(z dependence), this sample is combined with two 
other complete radio-selected samples which together 
span a wide range in L T &d at a given redshift; we also 
make a comparison with the SSQs in an optically- 
selected sample. In Section 3 we present evidence for 
a radio-optical correlation, and assess the role of sam- 
ple selection. In Section 4 we attempt a physical in- 
terpretation of the radio-optical correlation, and com- 
ment briefly on the implications of this result for stud- 
ies of radio galaxies. Throughout this paper, we take 
values of Ho = 100 kms~ 1 Mpc _1 , Q, = 1 and a zero 
cosmological constant. 



2 MOLONGLO-APM QUASAR SURVEY 

Our new, complete sample of SSQs was drawn ini- 
tially from the 408-MHz Molonglo Reference Cata- 
logue (MRC, Large et al. 1981), down to a flux den- 
sity 6408 > 0.95 Jy with a radio spectral index cut- 
off q > 0.5 (where a — — dlog S v /d log v evalu- 
ated near 1 GHz). The MAQS survey area is lim- 
ited to a ~ 1-sr region in which both UK Schmidt 
APM (Automated Plate-measuring Machine) data 
and Texas Catalogue positions (Douglas et al. 1996) 
were available, bounded roughly by —35° 1z S ^ 
0°, 21 h < q < 5 h . Hence, the sample is called 
the Molonglo/APM Quasar Survey (MAQS; Serjeant 
1996, Maddox et al. in preparation, Serjeant et al. 
in preparation). Optical identifications were made 
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Figure 1: The radio luminosity (L408); redshift (z) plane for SSQs from the MAQS (squares), the MQS 
(stars), 3CRR (circles) and the BQS (crosses). The brightest 33% (excluding the BQS) have been plotted 
as open symbols. For 3CRR the 178MHz-750MHz spectral indices from Laing, Riley & Longair (1983) were 
used. The loci of sources with a = 0.85 and S 408 = 0.95 (MAQS/MQS limit) and 4.94 Jy (3CRR limit) arc 
shown as dashed lines. 



from the UK Schmidt bj plates, using APM classifica- 
tions to exclude "galaxy" identifications brighter than 
6j = 20.5. The quasar candidates therefore consisted 
of all "stellar" or "merger" identifications, and all 
"galaxy" identifications fainter than 6j = 20.5. Red- 
shifts have been measured (at the Anglo-Australian, 
William Herschel and Nordic Optical Telescopes) for 
all quasar candidates to the plate limit of bj ~ 22.5, 
yielding a complete spectroscopic catalogue of 159 
(steep- and flat-spectrum) confirmed quasars. Opti- 
cal spectra of all quasar candidates (including those 
rejected for lacking broad optical-UV lines) will ap- 
pear in a forthcoming paper. 

The deep plate limit is very important to this 
study: only then is the survey likely to detect essen- 
tially all the SSQs as far as z ~ 3.3 where the Lya 
line leaves the bj band. A recent spectroscopic study 
of a complete and radio-fainter sample of 7C radio 
sources (Willott et al. 1997) quantifies this by finding 
no SSQs optically fainter than the plate limit used 
for the MAQS (i.e. SSQs with bj > 22.5 comprise 
less than 5 per cent of the total population) , but that 
most are fainter than the B ~ 20 limit of studies 
based, for example, on POSS-I. 

We have combined data on our new sample with 



data on SSQs from the radio-selected 3CRR sample 
(Laing, Riley & Longair 1983). In addition, we aug- 
ment the MAQS sample with SSQs taken from the 
Molonglo Quasar Sample (MQS, Baker 1994), also se- 
lected from the MRC down to S408 = 0.95 Jy, and to 
the fej plate limit. Finally, we include the SSQs from 
the optically selected Bright Quasar Survey (hereafter 
BQS, Schmidt & Green 1983). Figure 1 shows 408- 
MHz radio luminosity as a function of redshift (L408, 
z) for the four samples; the MQS points are plotted 
in separate symbols where the identification is not 
shared by the MAQS. Taken together, these four sam- 
ples provide a wide dispersion in L408 at any redshift 
z, allowing us to distinguish cosmic evolution (z de- 
pendence) from radio luminosity dependence. 

We have split the combined sample into two flux 
density bins; open symbols are used for the 33% of 
sources with 6*408 > 3.23 Jy, and filled symbols for the 
remainder. This split was chosen to yield a roughly 
even fraction in the brighter bin throughout the red- 
shift range. Measured values of radio spectral index, a 
were used where available (100% of 3CRR and BQS, 
97% of MQS, 72% of MAQS), and a = 0.85 assumed 
otherwise. Some of the brighter SSQs, and nearly all 
the flat spectrum quasars, should appear in the Parkes 
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Figure 3: The variation of apparent magnitude with 
redshift for a quasar of fixed absolute magnitude. The 
two curves show the variation for a numerical integra- 
tion of the composite quasar spectrum, and the cor- 
responding variation with a power law optical quasar 
continuum (as discussed in the text). The power law 
model is the brighter at a redshift of z = 5. Both mod- 
els have been normalised to produce the same apparent 
magnitude at a redshift of 1.5. 

catalogue (Wright et al. 1991): for a source on the 
MAQS radio flux limit the Parkes catalogue should 
detect all the radio sources with a < 0.4. We are thus 
confident that flat-spectrum sources have been almost 
entirely excluded, despite our currently incomplete a 
information. 



3 RESULTS 

3.1 The radio— optical correlation 

In Figure 2, optical luminosity is shown as a function 
of z, again split into bright and faint radio subsam- 
ples. In order to calculate the absolute B-magnitude, 
Mb, we transformed 6j magnitudes to B magnitudes 
using B w 6 j — 0.14 and adopted typical optical spec- 
tral indices of 0.5 (the value derived e.g. from the com- 
posite 7C SSQ spectrum by Willott et al. 1997). A 
more careful treatment of K-corrections was not pos- 
sible since a significant fraction of the objects are from 
the literature, having redshifts but not spectropho- 
tometry. Typical uncertainties in apparent magni- 
tudes are ~ ±0.2 mag; similar contributions to the 
absolute photometric errors are made by the uncer- 
tainties in the K-corrections. (Error bars are omitted 
from the plots for clarity.) 

Two possible further sources of error are galactic 
reddening and emission line contributions to the bj 
flux (intrinsic reddening will be discussed separately 
below). The galactic reddening is typically Ab ^ 0.2 
magnitudes, and the variations in galactic reddening 
are expected to cause an additional relative photo- 



metric error of ~ 0.1 mags (de Vaucoleurs & Buta 
1983a, b). This variation is smaller than our typical 
photometric errors, so we do not correct for it here. 
The MAQS extends to galactic latitudes only as low 
as | & | = 30°, so any galactic reddening is likely to be 
more significant for 3CRR; this would only increase 
the significance of our results discussed below. Unfor- 
tunately not all the quasars from the literature have 
published line fluxes, so it is not possible to apply em- 
pirical emission line corrections to the foj fluxes uni- 
formly over the sample. However, it is well-established 
in radio quiet quasar surveys that emission lines make 
a small contribution to the overall photometric errors, 
e.g. Schmidt & Green (1983). Since the bj passband 
is probably wider than the B selection passband of 
Schmidt & Green (1983), we might expect the line 
effects to be even smaller. This is verified in Figure 3, 
where we have integrated the composite quasar spec- 
trum of Francis et al. (1991), appropriately redshifted, 
over the fej passband. Also plotted is the prediction 
for an a op t = 0.5 optical power law, which is the con- 
tinuum slope determined in the composite SSQ spec- 
trum of Willott et al. (1997). The deviations from 
the power-law model in Figure 3 are dominated by 
the slightly atypical (for SSQs) optical spectral index 
of the Francis et al. (1991) spectrum, except at red- 
shifts z > 3.3 where Lya exits the bj passband. Note 
the smoothness of the composite quasar curve. 

The loci of SSQs at three apparent magnitudes 
are shown as dashed lines in Figure 2; the upper line 
represents the magnitude limit of the BQS (B = 16), 
the middle line, the 6j = 20.5 limit above which 
optically-extended objects were removed from the 
MAQS, and the lower line the 6j ~ 22.5 plate limit of 
the MAQS and MQS. From Figure 2, any incomplete- 
ness in the MAQS caused by rejecting foj < 20.5 opti- 
cally extended objects is only expected to be present 
at z < 0.4. Broad-line radio galaxies from 3CRR are 
also optically extended on POSS plates, and have not 
been included. A point-like reddened 3CRR quasar, 
3C22, is plotted in Figure 2 as two asterisks joined by 
a line, the lower symbol representing the measured 
Mb and the upper representing the magnitude cor- 
rected for reddening (see Rawlings et al. 1995). 

A clear tendency is seen for the radio-bright 
quasars at any z to also be brighter optically. This 
trend cannot be explained solely on the basis of selec- 
tion effects since, for example, the 3CRR SSQs have 
no optical magnitude limit yet are brighter on the 
whole than the MAQS /MQS quasars. 

Comparing on Figure 2 the SSQs of MAQS/MQS 
with those of 3CRR (restricting both to 0.4 < z < 2 
to ensure well matched comparisons) , we find a mean 
apparent magnitude of 18.36 ± 0.20 for sources above 
3.23 Jy (open symbols), but 19.22 ±0.14 for the re- 
mainder; the null hypothesis of identical distribu- 
tions is rejected at ~ 99% confidence by both the 
Kolmogorov-Smirnov and Mann- Whitney tests. For 
the combined samples in the restricted redshift range 
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Figure 2: The optical luminosity (Mb), redshift (z) plane (symbols as Figure 1). The upper dashed line 
marks the BQS magnitude limit, the lower line the MAQS limit, and the central line indicates the approx- 
imate threshold of reliable star-galaxy separation. The reddened quasar 3C22 is marked as two asterisks 
joined by a vertical line (see text). 



0.4 < z < 2, we find the optical and radio fluxes are 
significantly correlated 99.9% confidence using 
Spearman's p); the same is true of the whole sam- 
ple (i.e. unrestricted z, Figure 4). This is the first 
evidence that the apparent radio-optical correlation 
is not dominated by a redshift effect. 

Figure 5 shows the L408-MB correlation implicit 
in Figure 2. The correlation in Figure 5 is significant 
at ^> 99.9% confidence (using Spearman's p), with a 
best-fit slope of —2.5 x (0.6 ±0.1) (in the log 10 -L408 
vs. Mb plane) and dispersion of ~ 1.6 optical mag- 
nitudes (using 3CRR and MAQS/MQS). However, 
this slope and dispersion may be biased estimators of 
the true values, because the density of points on the 
plane depends on the weighting imposed by the bi- 
variate (radio-optical) luminosity function (LF), and 
the sampled comoving volume. 

Thus, by comparing the 3CRR and MAQS/MQS 
samples we find evidence for the radio-optical cor- 
relation which (for the first time) is independent of 
cosmic evolution. Is the correlation present in any of 
the samples separately? It might be argued that if a 
correlation only appears when samples are combined, 
then it is more suspect than a correlation from a sam- 
ple in which selection effects are uniform across the 



sources and where fewer corrections need to be made 
to inter-compare the sources. In fact the correlation 
also present in the MAQS in isolation (to a similar sig- 
nificance level), but it is perhaps instructive to show 
circumstances in which this type of argument can be 
misleading. Suppose we had two samples of quasars: 
the first has z = 1 ± 0.1 and B = 18 ± 0.1; the second 
has z = 1T0.1 and B = 22T0.1. It would probably be 
impossible to demonstrate a radio-optical correlation 
with either sample in isolation, but once combined 
the correlation would be far easier to detect. 

The L408-MB correlation in Figure 5 is supported 
by the radio properties of the optically faintest SSQs 
(the band between the lower two dotted lines in Fig- 
ure 2). With the exception of reddened quasars like 
3C22, which appears well above this band once cor- 
rected for reddening (Rawlings et al. 1995), this re- 
gion includes no 3CRR SSQs. It is, however, well 
populated with MAQS/MQS SSQs. Although red- 
dened examples certainly exist (Baker & Hunstead 
1995, Willott et al. 1997), SSQs are typically red- 
dened much less severely than 3C22 (or are reddened 
so severely, Av S> 1, that they will be classified as 
radio galaxies even on the basis of K-band images), 
since 3C22 has no clear broad emission lines in the 
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Figure 4: The flux-flux radio— optical relation for SSQs (symbols as Figure 1). Note the horizontal offset 
between the filled (lower) and open (upper) symbols. 



observer-frame optical. This argues strongly against 
L4og-dependent reddening as a cause of the radio- 
optical relation, an explanation which would in any 
case require post-shock temperatures finely tuned to 
the destruction of dust in, and only in, the most lu- 
minous radio sources. 

This is not to say that reddening does not af- 
fect the correlation in any way; we only claim that 
reddening alone cannot cause it. For example, it re- 
mains possible that reddening may affect the scatter 
in the correlation (see Section 4.2), or that the degree 
of reddening may be linked to the optical luminosity 
e.g. via dust sublimation. 

Finally, although our samples have sufficient cov- 
erage of the radio-optical-redshift parameter space 
to demonstrate a radio-optical correlation, they are 
probably not large enough to address any possible 
evolution in this correlation. Interestingly, taking the 
higher redshifts in isolation the evidence for the cor- 
relation is more marginal. However, this may simply 
be due to lack of data, and the restricted dynamic 
range in radio power at high redshiftsj^j (see Figure 1). 
The same high redshift radio-optical behaviour is also 

* Such a selection effect may also be responsible for the 



seen in numerical simulations (discussed below) sup- 
porting this interpretation. Nevertheless, it remains 
possible that the correlation weakens or perhaps fails 
at high redshifts, implying a different mechanism for 
jet formation at high (z Z 2) and low (z ~ 1) redshifts 
(section 4). 

3.2 BQS Outliers 

The eleven BQS SSQs lie in an atypical region of the 
correlation. This apparent anomaly is also present 
in the radio-!/ [o in] plot (Rawlings 1994), so is un- 
likely to be (for instance) due to photometric er- 
rors in the BQS. There are two points to made 
about this feature, which we will support in the 
next section with numerical simulations. First, there 
is a selection effect which necessarily over-populates 
this region: the BQS, which covers half the sky, se- 
lects the rare quasars with the most extreme Mb 
at any z, favouring the high Mb side of a radio- 
optical relation with large intrinsic scatter. In con- 
trast, the deeper MAQS is limited to ~ 1 stera- 

apparent tightening of the Q p hot vs - Qbulk correlation seen 
at high-2 in 3CRR (e.g. Rawlings & Saunders 199f). 
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Figure 5: The radio-optical relation for SSQs (symbols as Figure 1). 



dian. Second, the BQS SSQs are typically ~ 10 
times brighter in the radio than the faintest end 
of the radio LF, Uos ~ 10 24 ' 5 W Hz -1 sr -1 , where 
the comoving space density is an order-of-magnitude 
higher (Dunlop & Peacock 1990). The lack of BQS 
SSQs with fainter radio luminosities is difficult to 
explain without either a radio-optical correlation, 
or a strongly luminosity-dependent quasar fraction 
(Lawrence 1991; Jackson & Rawlings 1997; Serjeant et 
al. in preparation). For the quasar fraction to explain 
the deficit of I0W-L408 BQS SSQs, the SSQ luminosity 
function would have to be non-monotonic (i.e. num- 
ber density must not be a strictly decreasing function 
of luminosity) . There is no evidence for this within the 
MAQS (Serjeant et al. in preparation), though it is 
possible that fainter radio samples may find such a lu- 
minosity cut-off (e.g. Willott et al. 1997). This leaves 
the radio-optical correlation as the more plausible ex- 
planation. Furthermore, the fact that four of the BQS 
SSQs are also 3C radio sources (although only two 
meet the more stringent selection criteria of 3CRR), 
and are thus among the most luminous radio sources 
in the sky, is clear evidence for the radio-optical cor- 
relation. Nevertheless, in section 3.3 we also show that 
the the top left hand corner of Figure 5 may be less 



well sampled for SSQs than other regions, so the evi- 
dence of a radio-optical correlation in this part of the 
radio-optical plane is less compelling; a much more 
convincing demonstration is found to follow from the 
lack of quasars in the bottom right hand corner. 



3.3 Simulated radio— optical relations 

As the previous discussion of the BQS outliers illus- 
trates, the passage of flux limits across the radio- 
optical plane make a qualitative understanding of the 
selection effects on the radio-optical relation rather 
difficult. To clarify the situation we made numeri- 
cal simulations of the data. Note that a quantitative 
comparison of the radio and optical properties of our 
samples would involve estimating the bivariate radio- 
optical luminosity function, which will be discussed 
elsewhere (Serjeant et al, in preparation); here we 
restrict ourselves to reproducing only the gross prop- 
erties of our samples. 

In Figure 6 we show the results of Monte-Carlo 
simulations of the MAQS, 3CRR and BQS samples. 
Points were sampled randomly from the three dimen- 
sional probability distribution defined by the bivariate 
(radio-optical) luminosity function. In Figure 6a, we 
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Figure 6a: Simulated data without a radio-optical correla- 
tion. 



assume there is no radio-optical correlation, i.e. the 
bivariate LF is 



$(1/408, M B , z) <x e(M B ,z)P(L 4 08,z) 



(1) 



with the optical LF, e, taken from Boyle et al. (1990), 
and the radio LF, f), taken from the "LDE" model of 
Dunlop & Peacock (1990). We renormalised the Boyle 
et al. LF to unity, i.e. 



1 



e(M B ,z)dM B = 1 



(2) 



Figure 6b: Simulated data with a radio-optical correlation 

at all z; thus we are only using the shape of the optical 
LF, and the cosmic evolution is determined by the 
radio LF alone. Note that the MAQS SSQs in the 
(Mb, z) plane of (a) follow the &j plate limit closely, 
whereas the 3CRR quasars have no optical flux limit. 
For simulation (b) (Figure 6b) we use 



$(1/408, M B , z) oc ,3(1/408, z)7(£408, Mes) 



(3) 



with the radio-optical correlation 7 modelled as 
Gaussian scatter of 1.5 optical magnitudes about the 
plotted full line, and again adopting the LDE LF f3 
from Dunlop & Peacock (1990). In both simulations 
we used a = 0.85 to convert from 2.7 GHz to 408 
MHz. 

In both (a) and (b) we sampled 35 points from 
3CRR, 118 from MAQS and 11 from BQS, integrating 
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$ throughout 24.5 < log 10 (P/(W Hz" 1 sr" 1 )) < 29, 
— 20 > Mb > —29. Simulation (b) reproduces our 
data qualitatively, including the BQS outliers; the 
quantitative differences (particularly in the BQS z 
distributions) are easily explained if the SSQ optical 
LF J $ dI/408 evolves differently to the total quasar 
LF (e.g. La Franca et al. 1994), or if the SSQ ra- 
dio LF differs from that of radio galaxies (Serjeant et 
al. in preparation), or perhaps to unknown selection 
effects in the BQS itself (Goldschmidt et al. 1992). 
The qualitative agreement was found to be robust to 
the assumed slope or dispersion in the radio-optical 
correlation. Simulation (a), however, is grossly incon- 
sistent with our data: a general problem for mod- 
els without a radio-optical correlation is their inabil- 
ity to explain why SSQs from a bright radio sample 
(e.g. 3C) are among the brightest optically, and (more 
marginally) vice versa (e.g. why the BQS contains 
high-L 4 08 SSQs). 

An interesting feature of simulation (b) is the lack 
of evidence for a radio-optical correlation at high red- 
shifts, despite the fact that a correlation is present. 
In the simulation this is due to the lack of dynamic 
range in radio power at high redshifts; the correlation 
has a wide intrinsic dispersion and is not detectable 
without a wide range in radio power to compensate. A 
similar behaviour is seen in our 3CRR/MAQS/MQS 
combined sample at high redshift (Figure 2), again 
where the radio power dynamic range is least (Fig- 
ure 1). Also, in both our data and in simulation (b), 
the 3CRR quasars lie slightly on the radio-bright side 
of the radio-opical correlation. This may be analagous 
to the BQS outliers (section 3.2): 3CRR selects the 
radio-brightest quasars at any z, favouring the high 
1/408 side of a broad radio-optical correlation. 

Finally, a very clear (perhaps the clearest) visual 
demonstration of our correlation is obtained by com- 
paring Figure 5 with the Mb — £408 plane of simu- 
lation (a) (Figure 6a, lower panel). The 3CRR SSQs 
have no optical flux limit, so in the absence of a corre- 
lation there is nothing to prevent them having optical 
luminosities at the faintest end of the SSQ optical LF. 
This is clearly the case in simulation (a). However, in 
our data (and in simulation (b)) we find the 3CRR 
SSQs have optical luminosities lying on roughly the 
same radio-optical relation as MAQS/MQS. This is 
obviously difficult to account for unless the radio and 
optical luminosities are related in some way. 

In other words, the lack of SSQs in the bottom 
right hand corner of Figure 5 is real; this area has 
been well-sampled for SSQs. This is clearly highly sig- 
nificant, which we can quantify as follows. In 3CRR 
there are 31 SSQs in the range 10 26 < L 40 8 < 10 27 ' 5 
W Hz -1 sr~\ so if we assign absolute magnitudes 
uniformly over the range plotted in Figure 5 we find 
the probability that the area is empty to be about 
1CF 4 . This probability would be much smaller still if 
we were to use a more realistic optical LF. 

In contrast, the top left corner may arguably 



be less well sampled. The deficit of optically bright 
MAQS quasars could easily be attributed to the op- 
tical luminosity function of SSQs with faint radio lu- 
minosities. This leaves the BQS as the only survey 
which might adequately sample this region. The BQS 
does indeed lack the radio-faint, optically-bright SSQs 
which would fill this region, in agreement with what 
we might expect from the radio-optical correlation. 
However, the lack of these BQS quasars could per- 
haps be related to the incompleteness worries in the 
BQS (Goldschmidt et al. 1992). Several authors (e.g., 
La Franca et al. 1994) have also noted that the total 
quasar population appears to have a much higher frac- 
tion of SSQs than at higher redshifts. One suggested 
explanation of this, which may also help explain the 
apparent incompleteness in the BQS, is that at such 
low redshifts the host galaxies may contribute non- 
negligibly to the total magnitudes. As a result, the 
BQS (by eye) stellar selection may bias the sample 
with respect to host galaxy properties at these red- 
shifts in a complicated manner. 



3.4 Comparison with previous results 

At this point is is worth reviewing previous studies 
in the light of our correlation, and contrasting the 
much more problematic selection effects in previous 
samples. A common problem is the inability to dis- 
tinguish evolution from luminosity dependence. For 
example, in section 1 we discussed the apparent cor- 
relation in 3CRR between bulk kinetic jet power Qbuik 
and the photoionising radiation power output Qphot 
as estimated from the narrow line luminosity Lnlr 
(Rawlings & Saunders 1991). Also, Qbuik and Qphot 
have similar orders of magnitude in 3CRR, again sug- 
gesting a link. However, we will show that the selec- 
tion criteria of 3CRR could cause both this and the 
Qbuik-Qphot correlation, if 3CRR it taken in isolation. 

The Qbuik depends strongly on the total radio lu- 
minosity, which in turn correlates strongly with red- 
shift in 3C. This secondary correlation could ulti- 
mately lead to spurious relationships within 3C. For 
example, suppose there is no intrinsic Qbuik-Qphot 
correlation (implying that the correlation in 3CRR is 
due to some selection effect). Also, make the reason- 
able assumption that the optical luminosity function 
of radio-loud quasars evolves roughly as strongly as 
its radio-quiet counterpart. Then the narrow-line lu- 
minosity Lnlr will correlate strongly with redshift in 
3CRR, because Lnlr is evolving, and so the Qbuik- 
Qphot correlation in 3CRR would be due entirely to 
their independent evolution. The apparent similar- 
ity of the orders of magnitude of Qbuik and Q p hot 
throughout 3CRR may also be due to this indepen- 
dent evolution, and their apparent similarity would 
not be preserved in samples of fainter radio luminos- 
ity. A similar critique can be made of the radio-quiet 
Qbuik-Qphot relation, though in this case the samples 
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are optically selected rather than radio flux limited.|j 
However, the presence of a radio-optical correlation 
in our samples demonstrates for the first time that 
the Rawlings & Saunders (1991) result is not due to 
these selection effects, confirming their interpretation, 
and allows us to predict that the Qbuik — Qphot cor- 
relation will be preserved in complete samples with 
fainter limiting radio flux density. 

If we assume our observed dispersion in the radio- 
optical correlation is close to the intrinsic value, then 
we can also resolve some of the previous contradictory 
claims on the existence of the correlation. 

Peacock, Miller & Longair (1986) ruled out the 
null hypothesis that the apparently bimodal distri- 
bution of radio-loud and radio-quiet quasars is in 
fact due to a universal quasar radio-optical correla- 
tion (section 1). Our data can also immediately rule 
out the alternative model suggested by Miller, Pea- 
cock & Mead (1990), that all radio- loud quasars have 
a minimum optical luminosity of Mb — — 23. Also, 
on examining their sample selection it becomes clear 
why this study failed to detect the correlation. The 
Miller et al. RLQ sample was wisely restricted to a 
narrow redshift range, which counters any differential 
evolution, but unfortunately was limited to only six 
SSQs which spanned less than an order of magnitude 
in both radio and optical luminosities. Given the very 
broad dispersion apparent in our correlation, it hardly 
surprising that they did not detect it. Interestingly, 
although the Peacock et al. study did not explicitly 
address the SSQ radio-optical correlation, they noted 
that several of the quasars in their deeper Parkes sub- 
samples have lower optical luminosities than those in 
brighter Parkes subsamples. The samples were how- 
ever based in part on compilations from the Veron & 
Veron (1983) catalogue, and no attempt was made to 
separate SSQs from flat-spectrum quasars. 

Other previous inhomogeneous compilations also 
gave evidence for an SSQ radio-optical correlation. 
Neff, Hutchings & Gower (1989) found a clear radio- 
optical correlation distinct from evolution effects, 
though their sample was selected from the literature 
to fill the 2.7GHz radio luminosity, redshift plane 
as evenly as possible in 1 < z < 2. While this avoids 
the tendency towards fainter fluxes inherent in flux- 
limited samples, it is not clear if this method intro- 
duces selection effects of its own since the parent 
sample is clearly inhomogenous. Browne & Murphy 
1987 also found a clear radio-optical correlation in 
lobe-dominated quasars, though their sample was the 
radio-selected quasars in the Veron & Veron (1983) 
catalogue with published Einstein X-ray observations, 
which the authors emphasised "is a very heterogenous 
sample with all sorts of unknown selection effects." 

On the other hand, Browne & Wright (1985) had 



* Note though that the differing Qbulk/Qphot ratios in 
radio-quiet and radio-loud quasars is robust. 



the benefit of complete radio flux limited samples, but 
lacked complete spectroscopic redshifts. The radio- 
optical flux-flux correlation present in their Figure 1 
could then easily be explained by differing redshift 
distributions, instead of an intrinsic radio-optical cor- 
relation. For instance, one might reasonably expect 
that fainter samples extend to higher redshifts, so 
would have correspondingly fainter optical identifica- 
tions. Such an interpretation is ruled out explicitly in 
section 3.1 above. 

Miller et al. (1993) presented a study of the radio 
and optical properties of z < 0.5 BQS quasars, though 
the number SSQs was again too small to detect our 
correlation (see also Figure 5). However, we have also 
already noted that the completeness of the BQS has 
been questioned (Goldschmidt et al. 1992). 

In summary, the only previous SSQ samples with 
redshift- independent evidence for a radio-optical cor- 
relation were inhomogeneous compilations, so prone 
to unquantified (and probably unquantifiable) selec- 
tion effects. The complete SSQ samples, on the other 
hand, were either too small or too limited in radio 
or optical dynamic range to detect the correlation we 
have found. 



4 DISCUSSION 

4.1 A link between accretion and the 
fuelling of relativistic jets 

We argue here that the SSQ radio-optical correlation 
(as predicted e.g. in the jet formation models of Fer- 
reira & Pelletier 1995) suggests a close link between 
the formation of the jets and accretion onto the cen- 
tral black hole. Our discussion is similar to Rawlings 
& Saunders (1991), although now with the benefit of 
more direct evidence for a link between accretion and 
the fuelling of relativistic jets and without the selec- 
tion effect ambiguities. 

The narrow range in equivalent widths of broad 
emission features (e.g. Miller et al. 1993, Francis et 
al. 1993) and continuum variability studies imply that 
the bulk of the continuum is produced on sub-parsec 
scales, and is most naturally linked to accretion onto 
a black hole (e.g. Begelman, Blandford & Rees 1984). 
Therefore, the SSQ radio-optical link appears to arise 
on sub-parsec scales with the optical light produced 
by accretion, and with the radio luminosity derived 
from a centrally-formed jet. This confirms the view 
that the dominant influence on L408 is Q, and not 
the large scale radio source environment (Rawlings 
1993). 

It is of course possible that the accretion-jet link 
is not directly causal, since both processes could share 
a close link with a third parameter. Any correlation 
of the type shown in Figure 5 is often dismissed as a 
"brighter objects are brighter" effect. However, this 
scenario necessarily requires the existence of some 
links to create the scalings, and in the case of the 
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radio-optical correlation where both luminosities are 
generated by processes within the central parsec, such 
a link is very likely to be close. 

A radio-optical correlation could also be obtained 
if both Mb and L408 for individual quasars evolve 
separately, but in the same sense, with time. How- 
ever, if the radio lobes of SSQs are short lived (i.e. 
lifetimes <C Hq 1 , e.g. Alexander & Leahy 1987), re- 
quiring multiple generations of SSQs, then a conspir- 
acy would have to be preserved from 2 ~ 3 to 2 ~ 0.4 
despite strongly changing physical environments (e.g. 
Yee & Green 1987; Ellingson & Green 1991; Haehnelt 
& Rees 1993). 

The radio-optical correlation supports models in 
which SSQ jets are fuelled primarily by accretion on 
to a black hole (e.g. Begelman, Blandford & Rees 
1984), perhaps via magnetically driven winds (e.g. 
Spruit 1996; Ferreira & Pelletier 1995). In general, 
if any other parameters dominate the jet mechanism 
(such as disk angular momentum, disk structure, or 
related magnetic fields), then our correlation implies 
that they must also control or be controlled by the 
accretion rate. In particular, this allows us to exclude 
many classes of models in which SSQ jets are fuelled 
primarily by black hole spin energy and not accre- 
tion energy (Rawlings & Saunders 1990; Begelman, 
Blandford & Rees 1984). Only if the black hole spin 
regulates the accretion flow on (possibly) up to kpc 
scales can such models be sustained; Begelman (1985) 
presents a model in which he argues that local regu- 
lation of the accretion rate may be present. 

It is also worth remembering that radio-quiet 
quasars follow their own, but clearly different, 
radio-optical correlation (although it has not yet 
been shown to be redshift independent), having far 
lower jet powers than comparable radio-loud quasars 
(Miller et al. 1993; Lonsdale et al. 1995). Therefore, at 
least one of the above parameters differ in radio-quiet 
and radio- loud quasars. 

4.2 Scatter in the radio— optical correlation 

There are many possible physical interpretations for 
the broad scatter in our correlation. We will list some 
of the more obvious candidates here. 

First, radio luminosity gives an imprecise mea- 
sure of the bulk power in the jets (Qbulk) emanating 
from the central engine since it also depends on the 
gaseous environment of the radio source (Rawlings & 
Saunders 1991; Ellingson et al. 1991); a plot of Q vs. 
Mb may give considerably less scatter than Figure 5. 
This would be consistent with the smaller dispersion 
in the Q-Lnlr relation (Rawlings & Saunders 1991) 
for 3CRR radio sources. Unfortunately, at present we 
lack the radio data necessary to make the conver- 
sion from L408 to Q for the MRC and 3CRR SSQs; 
the compact, steep-spectrum minority may move the 
most in the conversion to Q because their strong con- 
finement and/or young ages make the L40S-Q con- 



version extremely important. These quasars may also 
obey a different radio-optical relation if they have sig- 
nificantly fainter or redder continuua than the main 
population (e.g. Baker & Hunstead 1995). 

Second, reddening of the optical quasar light is 
very likely contributing to the scatter (Baker & Hun- 
stead 1995; Baker 1996). It may be possible to reduce 
this cause of scatter by deducing reddening from the 
optical spectra, but since this typically needs higher- 
quality spectra and/or near- infrared photometry for 
our MAQS sample we have not yet been able to in- 
vestigate this. 

Third, optical variability may also increase the 
observed scatter and may also depend on luminos- 
ity. However, the variations are expected to be small 
for SSQs, about 0.2 mag on average on timescales of 
years (Hook et al. 1994). Variability scatter would be 
reduced by using median magnitudes over a longer 
timescale. The intrinsic spread in spectral energy dis- 
tributions over the observed optical waveband must 
also contribute to the scatter (e.g. Elvis 1994). 

Fourth, the sampling of the radio-optical plane is 
strongly non-uniform, as discussed above, which may 
be responsible for some of the apparent dispersion. 

Last, the dispersion could reflect additional de- 
pendences on secondary parameters. If this final sug- 
gestion is correct, then other observable quantities in 
SSQs may correlate with their deviation from the best 
fit radio-optical correlation. Such a discovery may 
provide far stronger observational constraints on the 
formation of the jets. 

4.3 The optical properties of radio galaxies 

The radio-optical correlation we have established for 
SSQs has important implications for the orientation- 
based unified schemes for active galaxies (e.g. Barthel 
1989; Antonucci 1993). The radio-optical relation 
suggests that amongst objects with a naked quasar 
nucleus those with the most luminous radio sources 
are necessarily also highly luminous in the optical. If 
the compact quasar nucleus is hidden from our di- 
rect view, as is now known to be the case in some 
radio galaxies (e.g. Dey & Spinrad 1996), then in the 
context of the orientation-based unified schemes the 
radio-optical relation predicts that a high narrow-line 
luminosity Lnlr is inescapable. The radio-optical re- 
lation is then an obvious candidate for the cause 
of the relation between narrow-line luminosity Lnlr 
and L408 for radio sources (Baum & Heckman 1989; 
Rawlings & Saunders 1991; McCarthy 1993; Rawlings 
1993; Jackson & Rawlings 1997). 

The SSQ radio-optical relation may also under- 
pin many anomalous properties of the most luminous 
high- 2 radio galaxies. For example, most models for 
the alignment of the radio and optical structures of 
high-2 radio galaxies (e.g. McCarthy 1993) require 
that roughly constant fractions of Qbulk and Mb are 
used to supply the different aligned components. If 
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we adopt the view that Q and Mb are closely linked, 
then as L408 drops so must the level of aligned light, 
whether it is supplied by Qtuik from the jets (e.g. 
Lacy & Rawlings 1994) or by Mb from the quasar 
(e.g. by scattering; Tadhunter et al. 1992). The lack 
of alignments and the low fraction of scattered opti- 
cal/UV light in less-luminous radio sources (Dunlop 
& Peacock 1993; Cimatti & di Serego Alighieri 1995), 
independent of redshift, support this prediction. 



5 CONCLUSIONS 

Using complete samples of SSQs, we have shown that 
the radio and optical luminosities of SSQs correlate 
with 3> 99% significance, independent of survey se- 
lection and cosmic epoch. Such a correlation would 
naturally explain the observations of aligned optical 
emission in only the most luminous radio galaxies. 
Also, we infer that the radio jets of SSQs are regu- 
lated by at least one parameter which is shared with 
the production of the optical continuum; in the ac- 
cepted standard model for active galactic nuclei this 
implies a link between accretion and the fuelling of 
the relativistic jets. 

This is not the first observational suggestion of 
such a link, nor the first claim of a radio-optical corre- 
lation. However, this is the first evidence that neither 
the correlation nor the apparent radio-optical links 
are caused by selection effects, particularly those in- 
herent in single flux-limited samples such as 3CR. 
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